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1. Introduction 

An application of effective field theory methods, which have become commonplace in 
high-energy physics and nuclear physics, has been achieved very recently by Brown 
and Yaffe [ij to analyze some general equilibrium properties of non-relativistic classical 
plasmas. In an exhaustive work, these authors have systematically derived a couple 
of higher-order results concerning with the thermodynamics and have computed some 
correlation functions of the three-dimensional classical plasma. Although of some of 
those results were known, the computation is streamlined by the organization inherent 
in the technique of effective interactions. 

The two-dimensional electron gas with 1/r interactions is both a fundamental model 
in many-body physics and a system of current experimental interest. Electrons trapped 
on the liquid-helium surface or electrons confined in the vicinity of a junction between a 
semiconductor and insulators or between layers of different semiconductors are examples 
of two-dimensional electron systems [2]. At finite temperature, the parameter nX^ 
measures the dilution of the system and the dimensionless coupling parameter defined 
hj g = 27™/5^e^ characterizes the strength of Coulomb interactions, where A denotes 
the thermal wavelength and e is the effective electronic charge in unrationalized units. 
For electrons on a helium film the electron density can be varied in the range from 10^ 
cm~^ to 10^ cm~^. For these low densities, it appears that the electrons form classical 
two-dimensional systems even at temperatures of a few millidegrees Kelvin. For T ~ 1 
K and 10^ < n < 10^ cm~^, the range of the plasma parameter is 1.7 < g < 17, not 
very far from the domain of weak coupling. 

In this paper, we largely pursue the methods advocated by Brown and Yaffe, 
and consider a classical two-dimensional multi-component plasma interacting via the 
Coulomb potential. From the appropriate scalar field theory, we compute the partition 
function and the correlation functions to first nontrivial order requiring renormalization. 
The analytical results that we have derived are valid in the low density and weak 
coupling domain where nA^ -C 1 and -C 1. The essential ingredients in the calculation 
are the matching coefficients or induced couplings, whose finite parts are determined 
from a comparison of the short-distance behaviours in the effective theory and in the 
full quantum theory. The main result of this paper is precisely the computation of 
the matching coefficients, exploiting some properties in the momentum space of the 
Coulomb Green's function for the quantum problem. 

There are old analytical and numerical results on the two-dimensional one- 
component plasma in the classical regime |Sl IH El Ej ■ Chalupa jlj reports an equation 
of state given by 



and Totsuji |5j gives the ratio of the correlation energy density to the kinetic energy as 



(3p 



I In + I (27E -l + \n2) + 0{g^ In" g,gHn g, g') , 



(1) 



n 



f3E, f3u 



l = |ln(7+| (27E-l + ln2) + .... 



(2) 



n 



n 
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These expressions agree precisely with our results (^Hjl and (pn|) for the pressure and 
the internal energy when the ratio of the Coulomb energy of a pair of charged particles 
separated by a thermal length to the thermal energy is very large, rj = f3e^/\ ^ 1. 
However, as we will show below, our one-loop results do not imply restriction to this 
ratio. They are valid to order g for arbitrary rj. The terms of order g'^ retained by 
Chalupa and Totsuji are not reliable because they must be computed by matching with 
unknown results of a three-body calculation in the quantum theory. However, we have 
been able to compute all the (^"(In f?)" leading logarithmic terms of the partition function 
in (j37|) . 

The plan of this paper is arranged as follows. In section El we briefly review the 
structure of the effective theory, the partition function to one-loop order and write 
down the form of the induced couplings. In Section IHl with the matching procedure 
completed, we compute at one-loop level the equation of state, the internal energy and 
the number density correlators. The section |3] is devoted to report some higher order 
results concerning the leading logarithmic contribution to the partition function. A 
short conclusion is presented in section |S1 Details of the derivation of the basic formulae 
for the matching procedure are given in the appendix. 



2. One-loop divergences in the effective theory 

We consider a plasma of different species of charged particles interacting through the 
Coulomb interaction. The charge, mass,. . . of the species a are denoted by e^, ma,. . . . 
These are the conventions of P which we shall closely follow. In the classical limit, the 
free-particle distribution function reduces to a gaussian function = e^^^e"^^^/^*"", 

and the free-particle density in two dimensions is n° = gaXa'^e^^", where 

K^i'-^r (3) 

denotes the thermal wavelength of the species a and ga, fJ^a are the corresponding 
degeneracy and chemical potential. We shall assume that the system is electrically 
neutral. 

The Fourier transform of the 1/r potential in two dimensions is 

-ikr 2tt 

which gives rise to the interaction potential in v spatial dimensions 

As usual, in the absence of any scale, the analytically continued coincidence limit V^(0) 
will be taken as zero. Using the fact that y— V^/(27r) is the inverse operator of 1/r 
in two dimensions, the grand canonical partition function in the classical limit can be 
turned via a Hubbard- Stratonovich transformation into the functional integral 



27r 



|D0(r)exp(-5ei[0;/i]), (6) 
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with the action functional of the electrostatic potential (f) defined by 

= / d'r |^0(r)^^0(r) - E^a(^)e^^^"^(^)} • (7) 

where n°(r) = ga\a'^e^^'^'^\ The corresponding field equation is the analog of the Debye- 
Hiickel result in two dimensions. Here, it proves convenient to assume that the chemical 
potentials entering into ra^ can depend upon r. Thus, the functional derivation with 
respect to fiai'^) will produce the correlation function between densities 

Kabir - r') = "f' 8 

Since the total charge neutrality in terms of the free-particle densities is 

E = 0, (9) 

a 

the = configuration is a solution of the classical field equation and, consequently, 
the functional integral can be perturbatively computed from a saddle point expansion 
around this trivial solution. The quadratic part of the action. 



/3 



2 



27r 27T 



(10) 



includes the inverse Debye length, kq = 2TTi3J2a ^a'^a? ^'^d the remaining part giving rise 
to the perturbative expansion is 

A5[0; /i] = - / d'r n° {e^^-<^W - 1 + ^/32e^0(r)2} . (11) 

This effective field theory can be systematically obtained from the quantum 
statistical mechanics of the plasma by integration of the quantum fields with momentum 
scales of order or larger pi. In addition to the terms already present in (jZj), 
there are sub-leading contributions in a derivative expansion in powers of (A90). The 
first of these comes from the static non-interacting two-point charge density-charge 
density correlation at non- vanishing k and zero frequency, evaluated within the Maxwell- 
Boltzmann statistics. This reads 

n(fc,a; = 0) = ^fl-^ + ...V (12) 



247r 

where the piece independent upon k produces the Debye screening included in (|Zj) and 
the next term yields an action contribution linear in the particle densities, 

•^V0(r).V0(r). (13) 

This action must be made consistent with the invariance of the original theory under 
the simultaneous constant shift (p (p — ic and fJ^a fJ^a — ^aC- Consequently, the 
corresponding derivative interaction, linear in the particle densities must have the form 

Sirl[<P;t^]=T.J rf^r^[V(/i,(r) + ze„0(r))]^n°(r)e^'^-<^W, (14) 
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where we have used the form of the Debye wave number. 

In the following, we safely use the effective theory of equation (|7j) to describe 
distance scales of order A and larger. Although this classical regime is independent 
of the quantum statistics, the logarithmic ultraviolet divergences in the effective theory 
are related to quantum properties of the full theory. These divergences of the form 
ln(fi;o//i), where /i is a momentum scale introduced by dimensional regularization, must 
match onto counterterms proportional to ln(/iA) arising from the action with the induced 
couplings required to cancel the fj, dependence. The precise value of the finite part of 
the induced couplings is determined by comparison with the result of the computation 
of some quantity in the full quantum theory. The main result of this work is justly the 
explicit evaluation of the one-loop induced couplings. 

The dimensionless loop expansion parameter is the ratio of the Coulomb energy 
for two particles separated by a Debye length to their average kinetic energy in the 
plasma. This plasma parameter is g = f3e'^KQ. Apart from the Debye length, there is 
another relevant length scale, not directly entering into the effective theory, but playing 
an important role in the matching procedure. This is the Bohr radius {e'^mab)~^ or, 
equivalently, the binding energy e'^mahl'2. of two particles in the plasma with reduced 
mass niab and equal and opposite charge. 

The Green's function Gy{r — r') in u dimensions, including the effects of static 
screening, is given in terms of Struve functions H and Bessel functions of second kind 
Yhy 

d^k e^'=-(^-^')27r 



G.(r-r') = / 



(27r)^ k + Ko 



AV,i\r-r' 



with 



AVJr) 



cos(7rz//2) 



Hi_w2(/«o^) - y\-vi2{i^^r) 



and the coincidence limit is 



G.(0) 



-vl2^v-\ 



(15) 



(16) 



(17) 



r(z//2) sinvrz/ 

With these results in hand, the one-loop contribution to InZ/A for the effective theory 
is given by 



E 



1 

2A 

Attu 



InDet 



G,AO) 



1 



where the determinant has been evaluated by integration of the formula 



(5 InDet 



1 + 



1 



d'^rG^iO 



5ko 
2tt ■ 



^..„j j--^..^,,^- (19) 

This contribution is singular when z/ — 2 and, already to this order, renormalization 
is required. To determine the form of counterterms, we note that the pole term is 
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a b c 



Figure 1. One-loop diagrams contributing to the irreducible part Cab{k) of the 
density-density correlator. The diagram c represents the contribution of induced 
couplings. 



proportional to Kq. Thus, at this order, the induced couphngs must be quadratic in 
the particle densities. The counterterm action, to be handled at tree level to cancel the 
one-loop divergences, must have the form 

/^] = E / /5^^7»?e^'^-*We^^-^W, (20) 

a, 6 

where the pole term of (/^^ 

(21) 



9l = t^'-' 



+ 9abW 



2 2-u 



cancels the u —>■ 2 pole of In^ coming from 6*2(0). The momentum scale /i is required 
in order for the finite part of the coupling QabifJ') preserve its dimensions at u = 2 when 
u is arbitrary. To this order, the contribution of the induced couplings only involves 
the = part of S'^d but at higher orders, when it is neccesary to handle the induced 
interaction beyond the tree level, the full exponential dependence upon (p will be crucial. 
The renormalization group equation for gab{fJ^), 

^^il^^'^^^ = (22) 

guarantees that the bare induced couplings g^f, do not depend on /x. By integration, one 
obtain 

gM = ^\^(4-]^ (23) 



where the integration constant fiab cannot be determined within the effective theory. The 
simplest way to determine the short- distance parameter fiat is to match the quadratic 
part in the particle densities of the one- loop density-density correlator Kab{k) in the 
effective theory and in the full quantum theory pp. In this comparison of the short- 
distance behaviours, the Debye screening plays no role P| provided that the correlators 
are evaluated at non-zero wave number k. 

As discussed by Brown and Yaffe , the Fourier transform of the complete density- 
density correlator has the general form 



I3{\k\ +2tI (3 Y.mn(imenCmn{k)) 
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Here, Cab{k) denotes the Fourier transform of 



Cabir 



(25) 



0=0 



where r[0; /x] is the effective action of the theory. At tree level, r[0; /i] = Sci[(p; fi] and, 
consequently, C^,j^'^{k) = Satn^ (no sum over a). Noting that 5'j^j[0;/i] + Sl'^\[(f); fi] in 
equations (fT^ and must be handled at tree level, the contribution to Cab{k) up to 
one-loop order is given by 

^1 



^Df(fc)-2^: 



ab 



-PSabeXG.iO 



-2pXbT.9'bM-^S.bnl 



(26) 



where 



Dl^\k) = J rfVe-'*'-"G,(r)2. (27) 

The function D^\k) comes from the one-loop diagram a in figure Hand Gy{0) is the 
tadpole. Since the last term of equation (j^Bj) yields a contribution to the n^ji^ part of 
Kab{k) which vanishes linearly when k —>■ 0, this term can be ignored in the matching 
procedure at small momentum k. Therefore, ignoring Debye screening. 



D^^\k) Cj^\k) = / d-'re 



ikr 



>2--7ri-f r(i 



K(r) 



f)r(V)^ 



■k 



u-2 



i2t 



This replacement leads to the piece of the correlator to be compared to the corresponding 
result in the full quantum theory 



nr. 



p2 2 



-Ci^\k)-2gl, 



27rf3nln°,ea 



k 



7e + In 



+ 2gabi^) 



27: (inyilegeb 
k 



(29) 



3. Matching and one-loop complete results 

In order to perform the matching, we consider the — * behaviour of the second- 
order contribution in the fugacity expansion to the quantum mechanical density-density 
correlator }Cab{k). The details of the computation are given in the appendix. This 
contribution to the full correlator is 

/Ci?(fc) = nyXb [F4k) ± {6ab/9a)F4k)] , (30) 

where the ± sign accounts for Bose (+) or Fermi (— ) statistics, ga denotes the spin 
degeneracy of the species a, and Xab is the thermal wavelength for the reduced mass 
i^ab = ''^a^b/ (^a + "^b)- The required terms are (Fourier transformed) matrix elements 
of the heat kernel constructed in terms of the Hamiltonian H for relative motion, and 
read 

F±(fc) = J dVe-*'=-^(r|e"^^'"'| ± r), (31) 
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with 



H, 



ab 



+ 



2mab r 

With the results (TOol) and (TOsll in hand, the comparison of ()30|) with 
enables to determine the integration constant fiab in the running coupling 



(32) 



gabifj^) 



(33) 



[ln(/i^A^J + f{r]ab) ± iSab/9a)fiVao 

in terms of the quantum parameter rjab = PeaCb/Xab- Thereby, the complete contribution 
to In Z up to one-loop order in terms of the fugacities is 



lnZ(/.) =J2n:-^Y.«{Pe.eb) 



A 



In 



471 



- 1 + 7e + r 



ab 



, (34) 



where 

^ab = f{riab)±{,5ab/ga)f{Vaa)- (35) 

By inversion of the formula 
1 d lnZ(/i) 

the elimination of the chemical potentials in favour of the mean number densities Ua 
can be made perturbatively, and at one-loop order we obtain 



1 + - Jl'^f' (/^eaCfe)' 
^ b 



In 



Atx 



+ 7E + r 



ab 



(37) 



with TtQ = 27r/5 X)a Ca^a- The corresponding replacement in the expression for the 
pressure p = \n.Z/ (3 gives the equation of state 

TT 



a. ^ a,b 

The internal energy u density is given by 

d lnZ(/i) 



In 



47r 



+ 1 + 7e + r 



which yields 

a 

where 



u 



TT 



dp A 



- J2 ri'anb iPeaCb) 

a,b 



In 



47r 



- + 7e + r^fe + - ?7„br'„fe 



= f'ir^^b) ± (5ab/^?a)/'(r7. 



(38) 



(39) 



(40) 



(41) 



In the case of a plasma consisting of a single species immersed in a uniform 
neutralizing background, the terms of order and In e which appear in the equation 
of state reported by Chalupa |1] agree precisely with our result when the strong 
repulsion result ()A.39j) is inserted^ into (jHSj), neglecting the term of order r]'^ and the 
exchange term ()A.50|) . The ratio of the correlation energy density to the kinetic energy 
density (3Ec/n reported by Totsuji ^\ also agrees with the term of order in pn|) when 
^aa > 1 • 

I The parameter in T is 47rne^/3^ = ^K^Xaailaa, and the parameter e in (51 is 27rne^/3^. 
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The insertion of our findings into (j26j) yields the irreducible part of the number density 
correlator function to one-loop order, in terms of bare densities 



Using (fTTj) and (PT|) . the contact term reads 



J2 TT-b iPeaCb) 



In 



+ 7E + r 



ab 



247r ' 



(42) 



(43) 



and the counterpart 6abnaFa{k) in terms of the mean number densities is readily written 
with the aid of (jHTjl . 



Fa{k) 



In 



An 



+ 7e + r, 



ab 



■ (44) 



Since G2{r)'^ — V2(r)^ behaves as ln(/tor)/r as r - 
when 1/ — >■ 2 may be written as 

D^^\k) = Cl'\k) + ACi'\k), 

where AC2'\k) is the Fourier transform of 

G2{ry-V2{rf = AV2{r) i^AV2{r) 

and the asymptotic behavior of AV2(r) reads 



0, the non- vanishing part of D^\k) 

(45) 



AVoir) 



{Yo{KQr) - Ho(Kor)) 



— + 



2 /vi3 



+ 



With this definition, we obtain 

7.0 n„\ _ A ^(2) , 



i^aVfc) = ACr(fc)-7r 



In 



^abl^ 



647r 



+ 7E + r 



(46) 



(47) 



(4^ 



We must emphasize that the dependence upon In(fc^) is deceptive. In fact, an explicit 
calculation when fc — shows that 



ACf)(fc) = -2n 
and, consequently 



1 - In 



4«:o, 



Kq 



In 



k_ 



0{k' 



-TT 



+ 2 + 7e + r. 



ab 



(49) 



(50) 



in perfect agreement with the general result CabW = —dfia/dPfib- Note that, to one- 
loop order, the n^n^ piece of Cab can be written in terms of the mean densities by 
replacing n° — > n. 
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4. Renormalization at higher order 

The complete determination of two-loop and higher order contributions to In Z requires 
at least to perform a difficult three-body calculation in the quantum theory. However, it 
is straighforward to derive renormalization group equations for the leading logarithmic 
pieces of the renormalized multiloop induced couplings gij^,,,i^ that enter into the 
counterm action through 

Sinai'P; /^] = E E / d'rgl^.j'^^-'^nle^P^^^'^ . . . <e^^^^"^ (51) 

n=2 i\...in 

The coupling g^^,,,i^ is proportional to e^^"~^^ since at = 2, all factors of n° excepting 
one of them must be accompanied by a factor of /3^e^, in order to reproduce the right 
power counting. Thus, the induced coupling g'^^ . j^ first contributes to \q.Z/A at order 
n — 1 loop. This contribution is 

= EC.X---<- (52) 

-il...tn 

Another contribution at order (n — l)-loop which arises from the one- loop graph with 
one insertion of the induced interaction g,°, ,■ , is 

Itl-.-^"^""'" E ^°...._.<---<_.(e.+... + e.„_J^^G.(0) (53) 

i\...in-l 

- E 4....„-.< • • • + ■ ■ ■ + e.„_J^et 

h-.-in 

where we have only shown the pole part of Gy{Q). Since the bare coupling g^^ 
contains in general a pole in z/ — 2 of order n — 2, this term produces a pole contribution 
of order n — 1 to be cancelled by /("^i) . Moreover, there are other pole contributions to 

9ii...in ' ^ 

\nZ generated by {n — l)-loop diagrams without induced couplings and, consequently, 
the cancellation of pole terms requires that gf_^ have the general form 

gt-in = ^I^'^'^T;^ ^ (gl-in-Mh + ... + ei^_,fel + permutations) 

ly ^ ft 

" n-l 

(54) 



+ ^(n-l)(i.-2) 



E 1 2^k + gii-iniij) 



where i?,-'^'' comes enterely from contributions not included into /I" . Thus, the 
finite term at z/ = 2 of the Laurent series for the renormalization group condition 
dg^^ ^^/dn = gives rise to 

l^'^9ii...iM = -- [gh...in-i{eH + ■■■ + ei^-ife}^ + similar] + . . . (55) 

By defining the input value ga = —1/2, this equation recursively determines the 
leading logarithmic contribution to the induced couplings when fi is of order kq and 
the integration constants are of order A~^. Up to three- loop order, the leading-log 
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pieces of the partition function are 

a a,b 



A 



a,b,c 



\ea + 66 + Ccf + ... (56) 



a,b,c^d 

In the case of one-component plasma the series of leading logs precisely agrees 
with the asymptotic expansion for the incomplete r(0, z"^) function, with z — 
7rn°e'^/9^1n(KoA). The result is 



lnZ(/x) _ 

— ft 



oo 



A 



n / 



n 

y 



1 



l + -exp - r(o,- 



1, 



(57) 



that produces the leading-log pressure expressed in term of the physical number density 

(58) 



^ - 1 = 1^ _ ^2 ^ 5 3 _ ^ 137_5 _ ^ ^ 

where 'z — Trrie^/?^ In(KoA). Finally, within this approximation the internal energy 
density becomes 



n 



(59) 



5. Conclusion 



With Brown and Yaffc as a guide, we have shown how the technique of effective field 
theory can be used in order to compute the partition function of the classical two- 
dimensional plasma interacting via a 1/r potential. Renormalization in the effective 
theory is essential already to one-loop level, and two-body induced couplings are required 
to order e^. We have performed in detail the computation of the finite part of these 
couplings by exploiting the properties of the Coulomb Green's function in momentum 
space. We have derived explicit expressions for the one-loop equilibrium properties, 
without restriction on the quantum parameter r]f^ — PcaCb/Xab- The higher order 
induced couplings satisfy renormalization group equations which can be solved in the 
leading logarithmic approximation. In the case of one-component plasma, we have been 
able to derive the asymptotic series of leading logs for the partition function in terms 
of the incomplete Gamma function. 
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Appendix. Computation of F±{k) 

In the following computation of F±{k), we shall temporarily omit the indices a,b. We 
start with the evaluation of given by the contour integral of the Green's function 

G= (H-E)-^ 

F^ik) = [ l^e-^^ [ rfVe-'=-(r|-l-|r), (A.l) 
Jc 2m J H — E 

where the contour C encircles clockwise the cut along the positive real E axis and all 

the bound-state poles when < in the case of an attractive potential. To exploit the 

rotational symmetry of the problem and the properties on momentum space, it proves 

convenient to average over the angle of k and compute the spatial integral making use 

of 

/•oo \ 

/ drrjQ{kr)jQ{\p-p'\r) = -5{k-\p-p'\). (A.2) 
JO k 

This produces the expression 

dE 



where 



F+{k) = J^—e-^^G^{k,E), (A.3) 

GAk,E) = Jd'pd'p'{p\j^\p')^6{k-\p-p'\). (A.4) 

The key to the computation of G^{k,E) is the closed form formula for the 
momentum space Green's function G{p,p') = {p\{H — E)~^\p') of the two-dimensional 
Coulomb problem. This form has been given in following the clever treatment of 
Schwinger [7] . Let us summarize the main results concerning this Green's function. The 
momentum representation equation for the Green's function is 



(A.5) 



Based on the (conformal) correspondence between the points of the euclidean two- 
dimensional momentum space and the surface of the unit three-dimensional sphere 

Px = — o sin 6* COS (/), (A. 6) 

^ 2(cos^/2)' ^' ^ ' 

Pv = — : — , ,9 sin 6 sin 6, (A. 7) 

^ 2(cos^/2)' ^ ^ ' 

the element of area on the sphere and the two-dimensional measure in the plane are 
related by 

dn=(^^yd'p, (A.8) 

\Po+P J 



where Pq = y/—2mE is the positive real wavenumber corresponding to the negative 
energy E. The angular distance 7 between the orientations of the two points Q = {6, (p), 



Effective field theory of the classical two-dimensional plasma 13 

Q' = {9', (j)') on the unit sphere is related to the euchdean distance between the 
corresponding points p and p' by 

2sin;y = , , , , „A P-Pf- (A.9) 



Then, noting that the term \p — p"\^^ can be written as 

-^^=2p,{pi+p^r''\pl+p'''r''^Y.i^yr^^^ 

and the completeness relation of the spherical harmonics 

5{Vl - ^]') = (A.ll) 

l,m 

the Green's function is just 

G{p, p') = 8mpl {pi + p^nvl + P")-'" E ^^""fj^^^'^* , (A.12) 

l,m 1+1/2 

where 

v= A. 13 

is positive when the potential is repulsive. For an attractive potential (z/ < 0), the poles 
corresponding to the bound states are 

E, = (A. 14) 

2(/ + 1/2)2 V ; 

The next step is to express the remaining factor in the integrand of ()A.4|1 in terms 
of the spherical harmonics associated with p and p'. This is readily performed from the 
series of Legendre polynomials 

5(2 - 2COS7 - a^) = E —^Pi " yj ^/(cost), (A.15) 
with the identification 

where 7 is the angle between the orientations in the unit sphere corresponding to p and 
p'. Therefore, using the addition of spherical harmonics 

^Kcos7) = ^ E Yr{n)Yrm*, (a.it) 



one obtains the required formula 



2po 



2 



+ P 



6{Q-Q') + 0{k^). (A.18) 
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With these results in hand, we return to the computation of G+{k,E). Using the 
integral representation of Legendre polynomials ^U] 

Pi{cos-y) = — / (i0 (cos 7 + z sin 7 cos 0) , (A. 19) 

TT Jo 

and the hypergeometric series with b = u + 1/2 and z = cos 7 + i sin7 cos</) 



E 

1=0 



i + b 



F{b,l,b+l,z), {\z\ < 1) 



(A.20) 



the sum over I in the Green's function (jA.12j) can be easily written in an appropriate 
form to study the fc ^ or, equivalently, the 7 — > behavior. This form is based on 
the following transformation formula for F{b, 1,6+1, z), valid when | arg(l ~ z)\ < vr 
and |1 — 2;| < 1 (see, for example, p 559 of [lOj) 



Fib,l,b + l,z)=bY^ 



n=0 



[n + l)-V^(6+l)-ln(l-z)](l 



Therefore, by evaluating the (f) integral of the ra = term, we obtain 



P/(C0S7) 



-7e 



In 



7 



0(7ln7), 



(A.21) 



(A.22) 



where 7e is the Euler's constant and denotes the digamma function. Also, we need 

1 



^P;(cos7) = 
T^o 2sm7/2 



(A.23) 



Putting all pieces together, the non- vanishing portion of G+{k,E) as fc ^ is readily 
found by the integration of ()A.4|) . This yields, after subtraction of the = 
contribution, the result 



G+{k,E) 



2miy 
Pok 



mu 

Po 
2-ny 



27E-l + 21n( — 1 +2i)[- + v 



Xk{-t) 2{-ty 



27e - 1 + In 



.2 
4tt 



0{k) 



-ln(-t)+2Jl + ^^ 



0{k) 



(A.24) 



where t = jSE, and the dimensionless parameter y is the ratio of the Coulomb energy 
for two particles separated by one thermal (reduced) wavelength to the temperature 



Ipme'^ 



^ Vab 



PeaCb 



(A.25) 



2n Xab 

Thus, the ± sign accounts for the repulsive (+) or attractive (— ) cases. Note that 
arg(— t) is purely real when t is on the negative part of the real t axis, so that, the 
convention is — tt < arg(— t) < tt when < a.Tg{t) < 27T. 

It remains to perform the required contour integral. The result Fj^{k) 
parametrically depends on and we will distinguish between the small and large \y\ 
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limits. For both the attractive and repulsive cases, the small \y\ << 1 limit can be 
obtained from the power series expansion of the ip function 



^(1/2 + = - 2 In 2 + £ CH (-1)^(2'^ - 1) 



„n— 1 



n=2 



, 1 

^l<2 



(A.26) 



provided that the contour C cuts the real t axis at the left oi t — — 47ry^, in the case of 
an attractive potential, or at the left of i = 0, in the case of a repulsive potential. The 
insertion term by term of this series in the Hankel's expression for r(z) 

-1-= [ A(_t)-e-*, (A.27) 
r{z) Jc27ri^ ' ' ^ ' 

and its derivative with respect to z 



r{z) 



produces the final result, after subtraction of the — contribution 



AF+(fc) 



2Try 7ry^ 



647r 



7e + In 



fiy) 



Xk 2 

where the function f{y) is defined by the power series expansion 

\n-i _ An'^y 



(A.28) 



(A.29) 



-2 + 



-iTiCW + 0{yY. (A.30) 



Note that here y can be positive or negative. 

Now, consider the asymptotic behavior when \y\ oo. For the attractive case, 
the leading contribution comes from the lowest bound state, as it is easily noted from 
the dependence on e~^^ . Hence, all we need is the residue — e^'^^^ of the integrand 
e~^G+{k, E) at the bound state pole t = —iTry"^. Thus, when y — >• — oo, 

AF+(A;) ~e^'^^', (A.31) 

and 

2 



m 



Try' 



(A.32) 



To obtain the leading asymptotic behavior in the case of strong repulsion, it 
proves convenient to study separately the contribution counterclockwise around the 
ramification point t = 0, and the contribution from the discontinuity of G^{k, E) at the 
cut along the positive real axis. As we shall see, the discontinuity yields exponentially 
small corrections. Using the identity 



Im-?/'(l/2 + i-u) = — tanh(7r-u), 
and noting that Disc \n.{—t) — —27ii, the discontinuity of G+{k,E) becomes 



(A.33) 



G+{k, E + ie) - G+{k, E - ie) 



iw'^y'^ 



-1 — tanh 



y 



2^,2 



2in^y 



+ 



iir'^y'^ 



Vi , 
1 — tanh 



, {t > 0)(A.34) 
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The last term in the second hne produces an integral well defined for y > 0, which can 
be studied by inserting the representation of tanhz as a series of exponentials. The 
resulting expression contributing to AF^{k) has the form 



--y'ti-irr^^e^p (-t- 

n=0 "'U t Y 



where, as y — »• cxd, the leading behavior corresponds to the first term in the sum, which 
has been evaluated by the steepest descent approximation. The first term in ()A.34jl 
yields 

- V r #^ = + - 7E - ln(5) + 0{5). (A.36) 



27r t2 s 

It remains to compute the contribution from the small circle around the origin. There, 
—t = 5e*^, with —71 < 6 < TT, and ln(— t) = ln6 + iO. Then, the contribution to AF+(fc) 
may be obtained by using the large z behavior of the ip function 

oo r> 

^(1/2 + ^ Inz + ^ (l - 2i-2") (I arg(z)| < vr) (A.37) 

n=l 

where are the Bernoulli numbers. The angular integral of e~*G^{k^E) produces 
singular terms which serve to cancel the linear and logarithmic divergences in ()A.36|1 . 
in addition to non- vanishing terms when 5 — > 0. Putting the pieces together, one 
immediately finds all terms of AF_|_(fc) non-exponentially suppressed in the large y 
limit, 



2 



47E - 3 + In —-^ 



^ (A.38) 



\k 24 

Note that the only contribution of the asymptotic expansion of the '?/'(l/2 + z) function 
comes from the logarithmic and the z"^ terms. With this result, the asymptotic behavior 
of fiy) as ?/ — > oo is given by 

f{y) ~ In (l6V) + 37E - 3 + (A.39) 

It remains to compute the exchange contribution. The nonvanishing portion at 
A; = of the exchange contribution is given by 

F-{0) = ^e-^^G_(0, E), (A.40) 

where 

G40,E)= Jd'pG{p,-p). (A.41) 

Noting that the orientations on the unit sphere of the points p and —p are Q = {6, 0) 
and Q' = {6,(j) + n) respectively, the angular distance between them is 26. Thus, using 
flA.17|) . the sum over m in ()A.12|) gives a factor Pi(cos2^). Now, the integration over p 
can be performed with the aid of the following integrals of Legendre polynomials 

r de sine Pi{cos2e)cos'^ 6/2 = (A.42) 
Jo 2/ + 1 
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and, finally, the resulting series is summed using the series representation 



E 



-1)' 



2Z + 1 + 2z/ 



4j 



1 V 

4+2 



Hence, subtracting the = contribution, (^-(O, E) reads 



G_(0,E) 



1 V 

4+2 



3 V 

4 + 2 



(A.43) 



(A.44) 



Since the exchange contribution involves particles of the same species, only the 
computation of -F_(0) in the repulsive case ?/ > is required. In this case, we can use 
the integral representation 



1 , e" 
du- 



2 Jo coshu/4 



(Re z>0) 



(A.45) 



When y -C 1, we insert the power series expansion of the exponential e and perform 
the u integral with the aid of 



du- 



2r(n+l)[C(n+l, l/4)-C(n+l,3/4)], (Ren > -1)(A.46) 



cosh u/A 

where ({s, a) denotes the generalized Riemann zeta function. Then, the evaluation of 
the contour integral yields the power series expansion 



F(0) 



' 2 



E 



:-l)"[C(l + n, 1/4) - C(l + n, 3/4)] f y^y' 



n+l 



(A.47) 



2r((3 + n)/2) V 2 ; 

Similarly to the function / for the direct part, it is useful to define the function f{y) as 



/(y)^-^F_(0) = i-4C 



TT^y 



0{y% 



(A.48) 



where C is the Catalan's number. 

Finally, the positive large y limit may be obtained by evaluating the u 
integral ()A.45jl using the Laplace's method of asymptotic analysis. Inserting the Taylor 
expansion of (coshM/4)^^ about m = 0, produces a convergent series in -u < 27r which 
can be integrated term-by-term. Therefore, the Hankel's formula leads to 

F_(0) ~ (A.49) 



and 



Ky) 



27r?/2 



(A.50) 



apart from exponentially small corrections. 
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